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ABSTRACT 

The BFS method for alloys is applied to the study of Co growth on Cu(l 1 1). The parameter- 
ization of the Co-Cu system is obtained from first-principles calculations, and tested against 
known experimental features for low coverage Co deposition on Cu(100) and Cu(lll). 
Atomistic simulations are performed to investigate the behavior of Co on Cu(l 1 1) as a function 
of coverage. 

INTRODUCTION 

Recent experimental analysis of Co growth on Cu(l ll) has raised interesting questions 
regarding the low coverage regime and the identification of the driving mechanisms that lead to 
the observed structures. The understanding of such phenomena would allow for controlling the 
structure and morphology of the growing films to achieve the desired properties. Some of the 
observed features include triangular-shaped islands with a height of two layers above the surface, 
decoration of the Cu step edges with mixed Co and Cu clusters, and monatomic-deep pools of 
vacancies in the surface [1]. We present preliminary results of a modeling effort based on the 
BFS method for alloys [2], meant to elucidate the main characteristics of the early growth stages 
of Co films on Cu( 111). Analytical calculations for the determination of the surface structure are 
supplemented with temperature-dependent Monte Carlo simulations. After testing the validity of 
the BFS parameterization (done with first-principles methods) by reproducing basic known fea- 
tures of the system, the methodology is applied to the study of multi-layer islands, their shape, 
composition and orientation. 

THE BFS METHOD 

The BFS method [2] is based on the concept that the energy of formation of a given atomic 
configuration is the superposition of the individual atomic contributions AH =Ze,. Each 
contribution e, is the sum of two terms: a strain energy, e,' s , computed in the actual lattice as if 
every neighbor of the atom i was of the same atomic species i, and a chemical energy, e, c , 
computed as if every neighbor of the atom i was in an equilibrium lattice site of a crystal of 
species /, but retaining its actual chemical identity. The computation of e, s , using Equivalent 
Crystal Theory [3], involves three pure element properties for atoms of species i: cohesive 
energy, lattice parameter and bulk modulus. The chemical energy, e, c , includes two BFS 
concentration-dependent perturbative parameters (Ac 0 cu and Ac u co) [2]. A reference chemical 
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Lattice 

parameter (A) 

Cohesive 
energy (eV) 

Bulk 

modulus (GPa) 

P 

a (A 1 ) 

/(A' 1 ) 

X (A) 

Co 

3.45803 

5.56374 

255.21 

6 

3.01887 

0.26166 

0.73527 

Cu 

3.63255 

3.55808 

141.79 

6 

2.87172 

0.27390 

0.76966 


Table I: LAPW [4] results for the lattice parameter, cohesive energy, and bulk modulus for the fee phases of Co 
and Cu, and the resulting Equivalent Crystal Theory (ECT) [3] parameters p, a, / and X (see text and Ref. 2). 

energy, e, c ", is also included to insure a complete decoupling of structural and chemical features. 
From the theoretical standpoint, the immiscibility of Co and Cu in the bulk presents a challenge, 
due to the lack of accurate experimental evidence that would help validate the parameterization 
used. However, although these binary phases do not exist, their properties can be readily 
computed via first-principles methods. Therefore, in this work, all the necessary parameters were 
determined using the linearized augmented plane wave method (LAPW) [4], by computing the 
equilibrium properties of Co and Cu in the fee phase, as well as the equilibrium properties of 
Co.iCu (LI 2 ), CoCu (Llo) and C 0 CU 3 (LI 2 ) ordered alloys. Moreover, this calculational scheme 
enables an accurate determination of the concentration-dependence of the BFS parameters A. 
The functions used in this work are given by Ac<>cu(*g>) = 0.045003 ln(jct«) - 0.047823 and 
AcuCo(4'c„) = -0.2101 1 3xc„ + 0.410351, where xa> (A'c,<) denotes the fraction of nearest-neighbors 
of a given Cu (Co) atom that are Co (Cu) atoms. Finally, the strain and chemical energies are 
linked with a coupling function g, which ensures the correct volume dependence of the BFS 
chemical energy contribution. Therefore, the contribution of atom / to the energy of formation of 
the system is given by 

e, = e, s + gfcf- £, c ") ( 1 ) 

Table I lists the necessary parameters for applying the BFS method to the Co-Cu system. 
The procedure used to determine the A(x) functions exceeds the scope of this paper and will be 
published elsewhere. We refer the reader to Ref. 2 for a detailed discussion of the BFS method, 
its definitions, operational equations and their implementation. 

RESULTS AND DISCUSSION 


Obtaining atomically sharp interfaces between different layers is one of the main problems 
in the preparation of thin films or superlattices. Intermixing is a common phenomenon, even 
between elements that are immiscible in the bulk. A typical example is shown in Fig. 1 , where 
the STM image shows a Cu(l 1 1 ) surface after growing 0.3 ML of Co at room temperature. A 
detailed study of the image reveals the presence of multilayer, mixed composition islands, with 
two or more atomic layers over the Cu(l 1 1) surface [5]. It is also found that there are pools of 
single-atomic-depth vacancies between the islands, absent in clean Cu(lll) surfaces. Other 
remarkable features include the fact that the total volume of the islands exceeds the amount of 
Co deposited, as well as evidence for the segregation of Cu to the island surface, as confirmed by 
ion scattering spectroscopy [ 6 ] and chemical titration [5] experiments. These features are 
highlighted in the inset in Fig.l, where an enlarged view of one of the islands clearly shows 
changes in the right side corners of the islands. STM experiments [7] also suggest that the 
islands include a subsurface Co layer. In addition, theoretical calculations using empirical 
interatomic potentials indicate that two-dimensional Co islands become unstable above a certain 
critical size and “explode” upwards, dragging Cu atoms from the surface layer forming mixed- 
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Cu(100) 

Cu(lll) 

Co(O) 

2.525 

3.356 

Co(S)Cu(0)i 

2.403 

3.176 

Co(S)Cu(0)r 

2.359 

3.081 

Co(lb) 

2.632 

3.360 

Co(2b) 

2.696 

3.415 


Table II: BFS total energy difference (in eV) 
between different configurations including a Co 
atom and a Cu(100) or a Cu( 111) substrate and 
the corresponding clean surfaces: (a) a Co atom 
in an overlayer site (Co(O)), a Co atom in a 
surface site with the ejected Cu atom in an 
overlayer site at (b) nearest-neighbor distance 
(Co(S)Cu(0)i) or (c) far from the Co(S) atom 
(Co(S)Cu(O),), and with the Co atom (d) one 
and (e) two planes below the surface. 



Figure 1: A 50x50 nm STM image of 0.3 ML Co deposited at 
room temperature on Cu( 111). The islands that appear on the sur- 
face are 4. 1 nm high over the Cu( 111) surface, corresponding to 
a double layer of Co. The inset shows an enlarged view of one of 
the islands, highlighting the fact that the composition of the 
islands is not homogeneous. The black hexagons correspond to 
vacancy pools in the Cu surface. 


composition islands several layers high [8]. 

From a modeling standpoint, the growth pattern of Co on Cu with increasing Co coverage is 
related to the ability of Co to interdiffuse in the Cu substrate. A simple atom-by-atom BFS 
calculation of several configurations including one Co atom deposited on a Cu(lll) substrate 
indicates that an exchange mechanism by which Co atoms substitute for Cu surface atoms is 
favored. The calculation also shows that a lower energy state is reached if the Cu atom is allowed 
to occupy an overlayer site away from the inserted Co atom (Co(S)). Table 2 lists the energy 
differences between a pure Cu( 100) or Cu(l 1 1) substrate and a state where (a) the Co atom stays 
in an overlayer site (Co(O)), (b) Co replaces a Cu surface atom which stays in a nearest-neighbor 
overlayer site (Co(S)Cu(0)i), (c) far from the inserted Co atom (Co(S)Cu(0)f), (d) in the first 
(Co(lb)) and (e) in the second (Co(2b)) layer below the surface. The noticeable energy gain in 
performing the exchange from Co(O) to Co(S)Cu(0) is a result of the lower surface energy of 
Cu than that of Co. For both surfaces, the Co(S)Cu(0)f is the configuration with lowest energy 
indicating that, in principle, Co atoms would not interdiffuse beyond the surface plane. A more 
detailed discussion concerning the different BFS energy contributions and their consequences in 
the deposition of Co on Cu will be presented elsewhere. 

Beyond the simple case of one adatom, Monte Carlo simulations provide useful information 
on the resulting behavior for higher coverages. In these simulations, dissimilar atoms are allowed 
to exchange sites at nearest-neighbor distance, as long as the Metropolis criterion is met, i.e., the 
exchange is accepted if it lowers the energy of the cell or it is accepted with a certain probability 
if it does not benefit the energy balance of the system. An indicative example is shown in Fig. 
2.a, which displays the final state of a Monte Carlo simulation. Unless otherwise noted, all 
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simulations where performed at room temperature in a rigid 
Cu lattice. In all cases it is observed that the ejected Cu 
atoms diffuse on the surface layer, forming clusters. The 
overlayer is denoted by a rectangle above the surface layer. 
Clearly, the difference in coordination between the (100) 
and (111) surfaces is small enough to ensure that the same 
behavior is to be expected for a Cu(l 1 1) surface, as shown 
in Fig. 2.b. 

With increasing Co coverage, the simulations show 
that Co atoms form patches in the surface and lb layers 
(Fig. 3.a-c). The formation of Co aggregates is expected, 
due to the immiscibility of Co and Cu. The interdiffusion to 
subsurface (lb) layers, however, deserves further study. 

These patches, with high surface energy, act as pinning 
sites for the Cu islands (Fig. 3.d), leading to the formation 
of bi-layer Co/Cu islands with one of the layers immersed 
in the substrate surface layer. It should be noted, however, 
that the process of Co diffusion on the surface and insertion in the surface layer is greatly 
enhanced with increasing temperature, as shown in Figs. 3.e-f, which repeat the results shown in 
Figs. 3.c-d for T = 1 K, where Co insertion in the surface layer competes with the formation of 
shapeless Co/Cu patches in the overlayer. Moreover, it is precisely the diffusion of Cu atoms 
along the edges of the forming Co/Cu islands that determines their triangular form, as seen in 
Fig. 4, where an initial state characterized by a Co/Cu island of undetermined shape evolves into 
nearly triangular islands. 

The small difference in energy between Co(S) and Co(lb) substitutions, indicated in Table 2 
for the case of one single Co atom, facilitates the interdiffusion of Co atoms in the substrate, 
unless alternative mechanisms are allowed that identify other energetically favored processes. 
Experimental results, usually available for high Co coverage, suggest the presence of bi-layer 
islands (i.e., two layers above the surface) which, naturally, should be preceded by the formation 
of single-layer islands. As long as atoms are allowed to populate only the first overlayer, the 
resulting islands adopt a triangular form where Co atoms cluster in the center of the island, 
surrounded by a Cu brim, as shown in Fig. 5. In general, the Co core also adopts a triangular 
form, inverted with respect to the overall island shape, thus maximizing the exposed Cu surface 
and therefore minimizing the surface energy of the island. 



Figure 2: Final state of 4 Co atoms 
(black spheres) in a Cu (a) (100) and 
(b) (111) surface, at T = 300 K. 



Figure 3: Final state of Monte Carlo simulations of a Cu( 111) surface with (a) 0.056, (b) 0.083, (c) 0.167, and (d) 
0.250 ML Co coverage. In every case, the initial state is given by similar cells where Co atoms are randomly 
dispersed on the Cu substrate. All simulations were performed at T = 300K. except (e) and (f) where T = IK. 
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Figure 4: Initial and final states of a Monte 
Carlo simulation at room temperature, 
corresponding to 0.25 ML Co coverage. 
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Figure 5: Final state of room temperature Monte Carlo 
simulations of a Cu( 111) surface with (a) 0.069, (b) 0.194, and 
(c) 0.313 ML Co coverage. In every case, the initial state is 
given by similar cells with Co islands. 


However, this is not the lowest energy configuration. A calculation showing the competition 
between single- and bi-layer (on the Cu(lll) surface layer) islands shows that there is a 
maximum island size beyond which the growth of a second layer is favored. Fig. 6 shows the 
BFS formation energies of single- and bi-layer islands with the same number of atoms. It is seen 
that while the BFS chemical energy favors the bi-layer island, regardless of size, there is a 
noticeable gain in strain energy that drives the transition from single- to bi-layer growth beyond 
a certain island size. The Monte Carlo procedure used in the preceding simulations does not 
allow for the growth of a second layer above the surface, as exchanges of pairs of atoms are 
considered one at a time. However, creating an initial state with some Co atoms on the second 
layer suffices to trigger the corresponding process. Allowing the Monte Carlo simulations to deal 
with the growth of a second layer on the Co/Cu islands leads to a similar mixed composition 
layer at the expense of the subsurface Co, thus leading to tri-layer Co/Cu islands where the first 
layer is immersed in the Cu substrate. As observed in the single layer case, the simulations show 
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Figure 6: BFS total, chemical and strain 
energies for different 1- and 2-layer islands 
with varying island size. Single- (bi-) layer 
islands are denoted with solid (dot-dashed) 
lines. 



Figure 7: Final state of room temperature simulations of a 
Cu(l 1 1 ) surface with (a) 0.299, (b) 0.396, and (c) 0.507 ML 
Co coverage. In every case, the initial state consists of Co 
islands. 


that each layer in the bi-layer islands is of mixed composi- 
tion, with Co in the center of the island, once again 
exposing as much Cu as possible (i.e., inverted Co 
triangles), as indicated in Fig. 7. At this point, it is clear 
that the critical size or, correspondingly, the critical Co:Cu 
ratio in the composition of the island can be reached if 
additional sources of Cu atoms - besides those replaced by 
Co insertions in the surface layer - are available. Fig. 8 
shows the initial and final states of a simulation where 
additional Cu atoms were initially randomly distributed on 
the Cu(l 1 1) surface, thus simulating the possibility of Cu 
atoms (decoupled from nearby vacancy pools, as seen in 
Fig. 1) diffusing along the surface. This additional amount 
of Cu clearly optimizes the energy balance of the resulting 
island, drastically inhibiting Co interdiffusion to 
subsurface layers, thus achieving a final shape and 
structure consistent with experiment. 

CONCLUSIONS 



Figure 8: Initial and final states of a 
room temperature simulation of a 
Cu(lll) surface with 0.17 ML Co 
coverage. 


The basic features of the process of island formation in the deposition of Co on Cu(l 1 1) 
were described within the framework of a simple modeling approach based on the BFS method 
for alloys: 1) intermixing of Co with the Cu(l 1 1) substrate, 2) formation of a Co surface layer 
on which mixed composition single- and bi-layer islands form, 3) evidence for a triangular island 
shape, 4) the existence of a critical island size, establishing the transition from a single- to a bi- 
layer growth pattern, 5) details on the distribution of Co and Cu atoms in the different layers, 
both in composition and location of the different atoms, and 6) hints on the temperature 
dependence of the above mentioned processes. The BFS analysis of the energetics of the system 
allows for the identification of competing driving mechanisms for the observed behavior. 
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